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ABSTRACT 

The mass distributions of dense cores in star-forming regions are measured to have a shape similar 
to the initial mass function of stars. This has been generally interpreted to mean that the constituent 
cores will form individual stars or stellar systems at a nearly constant star formation efficiency This 
article presents a series of numerical experiments evolving distributions of dense cores into stars to 
quantify the effects of stellar multiplicity, global core fragmentation, and a varying star formation 
efficiency. We find that the different evolutionary schemes have an overall small effect on the shape of 
the resultant distribution of stars. Our results imply that at the current level of observational accuracy 
the comparison between the mass functions of dense cores and stars alone is insufficient to discern 
between different evolutionary models. Observations over a wide range of mass scales including the 
high or low-mass tails of these distributions have the largest potential for discerning between different 
core evolutionary schemes. 

Subject headings: stars: formation — ISM: clouds — ISM: structure 



1. MOTIVATION 

The mass of a star is the single most important pa- 
rameter in determining how it will interact with its en- 
vironment, how long it will live, and the nature of its 
death. Therefore the distribution of masses of newly 
formed stars the initial mass function (IMF) — has far 
reaching implications for the evolution of the cosmos. 
The IMF is usually assumed to be universal, with a 
shape described by a powe r-law above ~ 1 M Q (see Scalo 
1986| |2005| |Kroupa||2002| ), and a log-normal below 
Chabrierfl2003[). 



sec 



Stars form from molecular clouds, and therefore their 
nature and characteristics define the initial conditions 
for star formation. The hierarchical density and ve- 
locity structure of molecular clouds is indicative of su 



frso nic turbulence (Larson| |1981 Falgarone & Phillips 
90| |Williams et al |2UUg| 7H owever, cold, dense regions 
quiescent gas, or "cores," found within Galactic clouds 
are beli eved to be the sites of future low-mass star for- 
mation flVryers et al.|[l983{ |Benson fc Myers||1989| |Ladd 
|et"aill99lj i. 

A large 1.3 mm continuum survey of dense cores in the 
pOphiuchus star- forming region by iMotte et all (|1998|) 
revealed a core mass distribution with a declining power- 
law slope similar to the IMF. This led the authors to con- 
clude that the cores observed in thermal dust emission 
are the direct progenitors of individual stars or stellar 
systems. The dense core mass function (DCMF) in other 
regions has also exhibited similarities to the IMF (e.g., 



Testi fc Sargent 1998; Johns tone et al.|2001||Reid fc Wil 



son||2006| |Enoch et al.||2007| |Nutter fc Ward-Thompson 
20071. Cores in the Pipe Nebula identified through dust 



extinction show a turnover at masses about a factor of 
3 higher than the IMF of the Trapezium cluster lead- 
ing to the interpretation that cores evolve w ith a nearly 
const ant star formation efficiency of ~ 30% ( Alves et al. 
20071. 



As first suggested by Vazquez-Semadeni (1994), the 



density probability distribution function (PDF) in 
isot hermal turbulent flows is e xpected to be log-normal 
(see Elmegreen & Scalo 2004 for a review) providing 
theoretical means to produce the low-mass end of the 
DCMF. The power-law tail of the DCMF can also be 
explained in terms of post-shock gas within a turbulent 

2002), or by deviations 
TM5| . In addition, the 



medium (Padoan fc Nordlund 



from isothermality (Scalo et al 



star formation efficiency of cores is theoretically expected 
to be nearly constant and lie between 30%-50% if out- 
flows from protostars are the primary mediating factor 
( |Matzner fc McKeepOOO] ). 

These results support a one-to-one or nearly one-to-one 
relationship between dense cores and the future stars to 
form from them. However, the similarity of the DCMF to 
the IMF remains the only piece of observational evidence 
for this kind of relationship. Meanwhile, there are several 
reasons to think that this one-to-one relationship may 
not hold. 

It is clear that some dense cores must fragment to pro- 
duce the large fraction of observed multiple stellar sys - 
tcms ( puquennoy fc Mayor|[l991| |Goodwin et al.||2007[ ), 
and it is possible that several fragments pe r core may be 
necessary to explain close binary systems (Sterzik et al 



2003 Goodwin & Kroupa 2005). Unfortunately, most 



observations of pre-stellar cores are limited to spatial res- 
olutions orders of magnitude greater than characteristic 
binary separations. 

It is also difficult to determine what fraction of cores 
identified in survey data will evolve into stars. In a re- 



cent study of cores in the Pipe Nebula, Lada et al. ( 2007 1 



y un- 
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found that the majority of cores are gravitational 
bound, with only the highest mass cores appearing des- 
tined to form stars. The most massive core in the Pipe 
Nebula, Barnard 59, is the sole active core in the nebula 
and harbor s an association of 
[eTaT1[2007] . 

Despite these open questions, the intriguing similarity 
between the DCMF and the IMF remains. Given the nu- 
merous ways in which a core could possibly evolve into a 



20 young stars (Brooke 
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young star or stars, we construct numerical simulations 
in an attempt to quantify the effects of different core 
evolution scenarios on the resultant stellar IMF. We fol- 
low this introduction with an outline of our methods and 
results in § [2] Section [3] discusses difficulties in compar- 
ing astronomical datasets and expands upon our results 
in consideration of observational constraints. We then 
conclude in § [4] and present a brief look toward future 
studies. 

2. SIMULATED EVOLUTION OF THE DCMF 

2.1. The Models 

There are several mechanisms by which dense cores 
may form from m ore diffuse molecula r gas including am- 
bipolar diffusi on (|Mousch oyias|[1 991[ ), thermal fragmen- 
tation (Jeans 1961| |Larson turbulent fragmenta- 
tion (Padoan ct al. 19971), or through triggering events 



(see, e.g., [Elmcgreen| |r9T 



Padoan fc Nordlund| ( |2002 1 
derive a functional form for the PDF of dense cores cre- 
ated by turbulent fragmentation that is in good agree- 
ment with observations. 

All our simulations begin with a DCMF generated from 
the PDF of their Equation 24. The role of turbule nce in 
creating de nse cores is still under question (see, e.g., Kirk 
et al. ||2007). However, the results of our experiments do 



not depend critically on the validity of turbulent frag- 
mentation since it is the difference between the DCMF 
and the IMF that we are testing here, not the correctness 
of any particular formulation of the DCMF. We can thus 
take advantage of the convenient analytical form for the 
DCMF provided by this formalism. 

The simulations begin with a DCMF having a log- 
normal peak at fi c — 1.3 M Q , a dispersion <j c = 0.37 dex, 
and a power-law tail at masses greater than a few Mq 
with an index x c — 1.3 (dN tx mT Xa dlogm ), nearly 



equal to the Salpeter slope of the IMF (1.35; Salpeter 
19551. This shape is achieved using an Alfven Mach 
number, Ma — 5, a kinetic temperature, Tk = 10 K, an 
average particle density, n = 1000 cm~ 3 , and (3 = 1.7, 
where Ef. oc is the turbulent energy spectrum. 

The models evolve cores with masses from 
log(M/M©) = —2.5 to 3.0 in increments of 
AlogM = 0.1 dex according to the prescriptions 
described below. Each core mass bin maps to a stellar 
probability distribution created numerically through 
1000 repetitions. The stellar probability distributions 
generated from each core mass bin multiplied by the 
corresponding value of the core PDF sum to create the 
final IMF. The characterizing parameters of the final 
stellar distributions are /i s , a s , and x s , analogous to the 
DCMF. 

One-to-One: This simple model serves as a comparison 
model for other evolutionary schemes and is labeled ref. 
Each core forms a single star with constant star forma- 
tion efficiency SFE — 0.3. This creates an IMF with 
the precise shape of the DCMF but shifted in log M by 
0.5 dex. 

Variable Star Formation Efficiency: This model called 
SFEVAR evolves each core into a single star with uni- 
formly random SFE G [0, 1]. 

Multiplicity: These models explore two different multi- 





TABLE 1 






Model Results 




Model 


X s (J s 


Ms 




(dex) 


(M@) 


Control Model 


REF 


1.3 0.37 


0.39 


Variable SFE Model 


SFEVAR . 


1.3 0.48 


0.63 


Multiplicity Models 


MULTl . . 


1.3 0.46 


0.29 


MULT2 . . 


1.3 0.40 


0.42 


Fragmentation Models 


FRAGPDF 


2.0 0.34 


0.30 


FRAGUNI 


1.3 0.51 


0.23 


Composite Model 


COMP . . . 


2.0 0.53 


0.34 



plicity scenarios motivated by observations. For model 
multI, each core is converted into a single, binary, triple 
or quadruple system with a probability of 57%, 37%, 4% 
and 1% respectively to produce stellar syste ms with mul- 
tiplicity in acco rd with solar- type field stars ( Duquennoy 
& Mayor 1991). The system components are assigned 
uniformly random mass ratios and then scaled such that 
the total mass in the system equals the total core mass 
multiplied by a constant SFE = 0.3. 

A mass dependency is embedded in model MULT2 
where only binaries are considered, and the probabil- 
ity for binarity increases linearly with mass from 10% 
for core masses < 0.03 M Q to 100% for core masses 
> 100 Mq. The mass ratio between pairs is uniformly 
random. 



Fragmentation: Two different fragmentation schemes are 
explored by these models. For model frag PDF, each core 
is fragmented into smaller cores repeatedly with masses 
drawn from the initial core mass PDF until no mass 
remains in the original core. This model mimics fur- 
ther turbulent fragmentation of observed cores. Model 
fraguni follows this same prescription, but draws the 
masses of fragments randomly with uniform probability 
between and the mass of the original core. 

Composite: The composite model, COMP, combines the 
evolutionary formulas of all the above models. It applies 
a random star formation efficiency, a multI multiplicity 
scheme, and fragmentation according to model FRAGPDF. 

2.2. Results 

Table [l] reports the results from all our models. All 
stellar distributions begin with a DCMF having x c = 1.3, 
a c = 0.37 dex, and /i c = 1.3 Mq. The slopes of the stellar 
distributions are derived from power-law fits at masses 
greater than /x s + 2a s , where the peaks and widths are 
measured directly. 

A variable star formation efficiency does not change 
the power-law slope of the resulting stellar distribution 
from the initial DCMF. The peak mass of the stellar 
distribution, fi s is lower than [i c by the expected factor 
of 2, and the spread in SFE broadens the distribution by 
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Fig. 1. — Distributions of original core masses (black dashed line) 
and resultant stellar masses for a simple one-to-one core to star 
relationship (bl ack), for a stellar distributio n with a multiplicity 
in accord with iDuquennoy & Mayorl (|1991k (red), and a binary 
probability increasing with core mass (blue). 



0.11 dex. 

Figure [l] shows the results of multiplicity models 
multI and mult2 in red and blue, respectively. Model 
REF is shown for reference. The effect of stellar multi- 
plicity has a negligible effect on the power-law portion 
of the IMF for either model. The creation of low-mass 
companions in model multI broadens the peak of the 
IMF by 0.09 dex and also shifts fi s toward lower mass by 
~ 25% in comparison to model ref. 

The redistribution of mass takes place mostly in the 
self-similar part of the DCMF for model MULT2. There- 
fore there is little change in the overall shape of the resul- 
tant IMF. The mean binary fraction for all stars created 
in model mult2 is ~ 30%. 

The results of our fragmentation models are shown 
graphically in Figure [2] Again, model ref is shown 
for reference. The cores of model fragpdf fragment 
into smaller cores with masses preferentially near fj, c . 
Therefore the higher mass cores have a higher number 
of fragments — up to ~ 100 for the most massive cores. 
This creates a narrower stellar distribution with a signifi- 
cantly higher peak compared to model REF and a steeper 
power-law slope by 0.7 in the index. 

The stellar distribution of model FRAGUNI has a power- 
law tail identical to the DCMF. The downward turn of 
the blue curve at the high-mass end in Figure [2] is due 
to the limited mass range over which we applied this 
fragmentation scheme. The randomly assigned fragment 
masses do not change the self-similar part of the stellar 
distribution, but the excess of low-mass stars generated 
by cores spanning the entire DCMF widens the resultant 
IMF by 0.14 dex. The mean number of fragments per 
core in this scheme is ~ 2. 

The composite model, COMP, combines stellar multi- 
plicity, core fragmentation, and a random star formation 
efficiency. While the preferred mass of the fragmentation 
creates a steep power-law tail as in model fragpdf, the 
broader peak due to multiplicity and a variable SFE cre- 
ate a shape consistent with the original core PDF over a 
wider mass range than fragpdf (see Figure [3] below) . 



3. DISCUSSION 

3.1. Observational Difficulties 

The different core evolution schemes produce stellar 
distributions that are clearly discernible in our theoreti- 
cal modeling. But these simulated data offer the luxury 
of complete control. When comparing real, astronomical 
datasets of cores and stars, several difficulties arise. 

3.1.1. Obtaining a Representative Sample of Cores 

The unbiased identification of an ensemble of dense 
cores within molecular clouds that will definitively form 
stars is a difficult task. Many automated core or clump 
finding algorithms exist that produce reasonable results 
from dust continuum, dust extinction or molecular lin e 
data (e.g. |Stutzki fc Guesten|1990| [Williams et al.|1994) . 
However, these methods otter no measure of systematic 
errors in identifying bona fide pre-stellar cores which may 
dominate the Poisson errors assumed in analyses. 

Once a core forms a protostar, it is not clear how the 
mass of the remaining core is relevant to the IMF. Cores 
with embedded stars can therefore be excluded from an 
ensemble with the use of sensitive infra red observations, 
e.g., using Spitzer (Evans et al. 2003). The remaining 
starless cores, however, may not all form stars in the 
future. 

Recent studies of the Perseus star-forming region show 
that starless cores tend to be less massive than cores with 
stars (although there exist ma ny low-mass cores that 
harbor embedded s ources) (see |j0rgensen et al. 2007 



Hatchell et al. 112007 1. This could mean that many of the 



low-mass cores included in their samples are transient 
structures or perhaps still accreting material. Molecu- 
lar line data can supplement dust maps to determine 
the gravitational boundedness of cores, and hence the 
likelihood of them eventuall y forming s tars. In a study 
of the Pipe nebula cores, Lada et al. (20071 find that 
only the ~ 25% most massive cor es are gravitationally 
bound, and Johnstone et al. (2000) find that a majority 
of cores in Ophiuchus are stable against gravity. How- 
ever, molecular line studies of NGC 1333 and Ophiuchus 
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Fig. 2. — Mass distributions from our fragmentation models. 
The black dashed and solid curves are the same as for Figure [l] 
The red distribution describes stellar masses resulting from core 
fragments with masses drawn from the original core PDF, while the 
blue distribution represents uniformly random fragment masses. 
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Fig. 3. — Comparison between simulated observations of the 
DCMF and the synthetic IMFs of several of our models. The error 
bars represent y/N measurement errors. 

find that cores in those regions extend ing down to masses 
of ~ 0.1 My are l ikely to be bound (|Walsh et al.||2007 
Andre et al.p007 l. 



Low-mass, starless cores that appear gravitationally 
unbound might also be bound by external pressure, and 
are therefore potentially stellar precursors. In this case, 
the possibility of different evolutionary timescales for dif- 
ferent cores in the sample may need to be considered (see 
Clark et al.|[2007| for details). 



3.1.2. Comparing DCMFs to the IMF 

It is impossible to measure a DCMF for an ensemble 
of cores as well as an IMF for the stars that formed from 
them. Variations in the DCMF {e.g., compare Motte et 
al. 1998 Nutter & Ward-Thompson 2007) and the IMF 



(IChabrier 2003) from region to region therefore must be 
considered when comparing these distributions. More- 
over, the total number of stars to form from a distribu- 
tion of dense cores can never be measured directly. This 
means that the height of, or total area under, the stellar 
IMF cannot be used to discern between core evolution 
models. 

Observations of cores in star forming regions are cur- 
rently limited to samples of less than a few hundre d {e.g., 
Enoch et al.|2007| [Nutter fc Ward-Thompson|2007l |Alves 
et al.|2007p . Using our core PDF, we generate simulated 
data for a sample of 300 cores. This number was cho- 
sen to give a number of stars per A log M = 0.2 dex bin 
near the peak of the DCMF comparable to modern-day 
observations. 

Figure [3] shows the comparison between these simu- 
lated observations and numerous results from our mod- 
els. The errors on the simulated data correspond to >/N 
from counting statistics. The synthetic IMFs are rescaled 
by a multiplicative factor and shifted by an additive fac- 
tor in logM that produces the lowest % 2 value between 
-0.1 and 1.1 in log (M/M ). All synthetic IMFs produce 
X 2 values less than one except model fragpdf for which 
X 2 = 3.6. 

Our models produce the largest variation in the param- 
eters ii s and a s , while only mass dependent core evolu- 
tion affects x s . Fragmentation affects the width of the 



IMF and tends to lower the characteristic mass of the 
IMF compared to a simple one-to-one core evolution sce- 
nario. Variation in the SFE broadens the IMF while the 
mean SFE for an ensemble of cores affects the charac- 
teristic mass of the IMF. Additionally, mass dependency 
may exist in either core fragmentation or in the SFE of 
cores. Therefore there is significant degeneracy between 
the nature of fragmentation and SFE in determining the 
shape of the IMF. 

To break this degeneracy, accurate measurements of 
DCMFs and the IMF over > 2 orders of magnitude in 
mass must be made, where the most leverage is achieved 
at the high and low-mass ends of these distributions. 
This data must then be supplemented with independent 
observations, such as comparisons between the total mass 
in cores and total mass in stars to constrain the SFE, or 
spatial correlation functions of stars and cores to con- 
strain the fragmentation. 

3.2. Observational Constraints 

The SFE of dense cores in Perseus has been recently 
deter mined to be between 10%-15% (J0rgensen et al. 
20071. This result implies that the horizontal shift be- 



tween measured DCMFs and the stellar IMF should be 
a factor of order 10. However, our models show that it 
is not the SFE alone that contributes to the horizontal 
shift between DCMFs and the IMF; fragmentation can 
also contribute. 

The SFE in model MULTl derived by direct comparison 
of the DCMF and IMF appears ~ 27% lower than the 
true SFE of 0.3 due merely to the minimal fragmentation 
required to reproduce the observed multiplicity of field 
stars in the Galaxy. For model fraguni the apparent 
SFE drops to 0.17. Fragmentation effects could therefore 
lead to an underestimation of the SFE by as much as 40% 
for a given measured DCMF. Alternatively, it is shown 
through model mult2 that the opposite systematic effect 
could be an issue if higher mass stars preferentially form 
multiple systems; i.e., an overestimation of the inferred 
SFE from a measured DCMF. 

Without considering the possibility of core fragmen- 
tation, the horizontal shifts between measured DCMFs 



and the IMF imply SFEs of > 5 0% in Ophiuchus flMotte 



et al 
Nebu 



|1998| |Andre et al . 
a^Alves et _al.||2007|7 
-Thompson 2007). I 



2007(1, to - 30% in the Pipe 
6% in Orion (INutter & 



to 



Ward-Thompson 2007). "This wide spread may be par- 
tially due to a varying degrees of core fragmentation in 
the different regions. However, it is not clear to what de- 
gree this discrepancy is due to observational biases. For 
instance, the poorer physical resolution of observations 
toward distant regions may artificially drive the DCMF 
peak to higher masses because of confusion, thus lowering 
the inferred SFE. 

The number of cores containing multiple embedded 
protostars within a star-forming region is related to 
the amount of core fragme ntation that has t aken place 
there. From the data of |j0rgensen et al. ( 2007 , Ta- 



ble 1), anywhere from 4% to 23% of cores have formed 
multiple protostars with separations between ~ 1500 AU 
and ~ 4000 AU. This is a lower limit to the amount of 
core fragmentation in Perseus primarily because ~ 40% 
of cores are likely to fragment on scales smaller than 
1500 AU to p roduce a multiplicity in acco rd observations 
of field stars ( Duquennoy fc Mayor|[l99T I . 
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However, for comparison, we calculate from our frag- 
mentation models an upper limit to the number of 
Perseus cores expected to be observed with multiple pro- 
tostars. Thirty percent of core s with masses between 
0.1 and 4.9 M ( |Kirk et al.||2006[ ) fragment into multiple 
cores according to the FRAGPDF model, while 42% of the 
cores from the fraguni model potentially produce cores 
with multiple protostars. 

4. CONCLUSIONS 

If the structure and conditions within molecular clouds 
indeed determine the stellar IMF, the similarity between 
measured DCMFs and the IMF is likely a reflection of 
this relationship. It is clear from past studies that some 
of the cores identified in star-forming regions are the di- 
rect progenitors of stars or stellar systems. However, 
our simulations have shown that the overall shape of the 
IMF is robust against different core evolution scenarios 
for core masses between ~ 0.5 and 10M©. In light of 
the uncertainties in identifying a representative sample 
of cores in astronomical data, a direct comparison be- 
tween DCMFs in this mass range and the IMF cannot 
alone imply any particular evolutionary path from an 
ensemble of cores to a future generation of stars. 

The peak mass in our simulated IMFs, fi s , is affected 
by both the SFE and the fragmentation scheme adopted. 
Estimations of the SFE by comparing the turnover in the 
DCMF to the IMF could be off by as much as - 40% 
without considering the possibility of core fragmentation 
(including multiplicity). The widths of our simulated 
IMFs, (7 S , are also affected both by variation in the SFE 
and fragmentation. Therefore, independent constraints 
on the SFE {e.g., through comparisons between total core 
mass and total stellar mass) or fragmentation {e.g., using 
spatial correlation functions of cores and stars) must be 
used to discern between these two effects given accurate 
measurements of DCMFs and the IMF over > 2 orders 



of magnitude in mass. 

Differences in the power-law tails of the DCMF and 
IMF are sensitive to mass dependencies in core evolu- 
tion. Using a fragmentation scheme with a preferential 
mass scale we achieve a power-law slope change of 0.7, 
though other fragmentation scenarios produce no notice- 
able difference. 



4.1. Future Studies 

The massive clouds identified through 8 ^m extinction 
against t he Galactic background — infrared dar k clouds 
(IRDCs; |Perault et al"1|1996[ |Egan et al-l|1998|) — offer a 



promising target list from which to derive DCMFs over 
a wide mass range extending to higher masses than have 
been previously measured. These sources, however, are 
typically distant making individual cores difficult to de- 
tect and resolve. ALMA may be needed to produce useful 
DCMFs from these sources. 

To constrain the low-mass end of the DCMF, it is nec- 
essary to make sensitive, high-resolution observations of 
the closest star-forming regions. Current facilities have 
the resolution to detect individual cores, and instruments 
such as SCUBA2 may extend considerably the DCMF to 
lower masses over wide fields. Together with molecular 
line data and complete surveys of late-type stars and 
brown dwarfs in a variety of star-forming regions, the 
statistics at the low-mass ends of DCMFs and IMFs may 
soon be sufficient to draw detailed conclusions regarding 
the relationship between the DCMF and the IMF. 
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